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We present results of a global tuning of general purpose Monte Carlo models to published measurements of 
IP — * jets at HERA and 77 — > jets at LEP and TRISTAN. The principle free parameters in the tuning are the 
simulation of the underlying event and the choice of photon structure. Several combinations of models are ruled 
out by the data. Some consequences of the tuned models at a future linear collider are discussed. 



1. INTRODUCTION 

A systematic examination of the manner in 
which data is described by Monte Carlo mod- 
els can give valuable insight into the underlying 
physics involved. However, many such compar- 
isons are to a single result, limiting the extent to 
which physics conclusions can be drawn and rais- 
ing the possibility that agreement is achieved in 
one cross section at the expense of another. The 
aim of this study is to tune, for the first time, 
general purpose Monte Carlo models to the ex- 
isting and expanding set of jet photoproduction 
data. To that end, we have compared model pre- 
dictions to published inclusive jet, di-jet and 3-jet 
data from the ZEUS @-[| and HI ||| experi- 
ments at HERA, the OPAL [p|JTT| experiment 
at LEP, and the TOPAZ || and AMY @ ex- 
periments at TRISTAN. The inclusion of both jp 
and 77 data (see figure |l|) is particularly signifi- 
cant for constraining the tuning. 

Our tuning currently focuses on the role of 
the so-called "underlying event" in hadronic jet 
production and the extent to which this can be 
described by perturbative QCD inspired models. 
We constrain the multiparton interaction (mini- 
jet) models contained within the Monte Carlos 
described below, and use these tunings to esti- 
mate backgrounds at future colliders. 

Such a tuning can also lead to a number of 
practical benefits. The method employed allows 
us to check the consistency of the data contained 
in the various publications, which ranges across 
different colliders, experiments, years, energies, 



kinematic regions etc. A better description of the 
data by the models can be expected to lead to a 
reduction in the systematic errors due to detector 
corrections in future measurements. 
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Figure 1. Illustration of the correspondence be- 
tween the different 77 and jp event classes. 



2. THE MODELS 

The models tuned were HERWIG 5.9 0], in- 
terfaced to the JIMMY eikonal model for 
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multiparton interactions, and PYTHIA 6.125 
Pg| , which contains a new (since version 5.7) 
mode aiming to simulate the virtuality of the 
bremsstrahlung photons ]l7). For HERWIG, 
a previous tuning to DIS data (CLMAX=5.5, 
PSPLT=0.65) flf] was used. In addition to 
the default mode, preliminary investigations were 
carried out on a modification to HERWIG, 
whereby the intrinsic transverse momentum dis- 
tribution of the partons in the photon takes the 
form of a power law rather than a gaussian. This 
modified version, referred to as HERWIG+fc f , 
was motivated by studies at HERA [jl9| and at 
LEP |^0| which indicate that it can lead to a bet- 
ter description of the data. 

The main free parameters investigated in the 
tuning are the description of the "underlying 
event" and the the choice of photon structure. 
Multiple hard scatters, illustrated in figure]^, can 
be enabled or not, and it is expected that some of 
the distributions studied will exhibit great sensi- 
tivity to this phenomenon. The other facet of the 
underlying event is the minimum transverse mo- 
mentum of the hard scatter(s), hereafter referred 
to as p™ m . Sensitivity to this parameter is greatly 
increased when multiple interactions are turned 
on. As one goes to lower pr, higher parton den- 
sities are being probed, leading to an enhanced 
probability that more than one hard scatter will 
occur. It should also be noted at this point that, 
whilst HERWIG employs a cluster fragmentation 
scheme, PYTHIA uses string hadronization. 

A number of different parameterisations of 
the photon structure function were investigated. 
Only leading order sets were used since the ma- 
trix elements in the Monte Carlo models are lead- 
ing order. The sets used were GRV 92 ]2l| , 
WHIT2 H, SaSld and SaS2d |§, and LAC1 
@. Throughout this work the GRV 94 LO Q 
proton pdf was used in the simulation of HERA 
data. 

Finally, the overall normalisation of the Monte 
Carlo was treated as a tunable parameter. This 
is justified within a range of around a factor of 
two or less because of the uncertainty in the scale 
of a s . 




Figure 2. An example of a multiple scattering. 



3. FITTING METHOD 

This work was carried out within the frame- 
work of a package developed by the HERA com- 
munity to permit easy comparison of data to 
Monte Carlo p6| . Extensions were made to fa- 
cilitate the inclusion of the 77 data. 

The procedure for finding the best fit, for a 
given set of parameters, of the Monte Carlo to 
the data was as follows. An overall % 2 per degree 
of freedom across all the distributions (some 50 
in all) was defined as: 

2 _ 1 <A (MC(i) - Data(i)) 2 



where Data(i) and MC(i) are the values of the dis- 
tribution in a given bin i for the data and Monte 
Carlo respectively. The sum runs over the total 
number of bins n in all the distributions. 

The aim was to minimise this x 2 > an d a good 
fit is when it is approximately one or less. The 
normalisation of the Monte Carlo was varied to 
find the best fit (The sum in equation (1) is di- 
vided by n — 1 rather than n to take into account 
the resulting loss of a degree of freedom) . For the 
HERA distributions, each data plot was allowed 
to float within the quoted correlated error (typ- 
ically of 15-20%) and the number of plots was 
then subtracted from n as well. Otherwise, all 
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systematic and statistical errors were treated as 
uncorrelated and added in quadrature. 

The distributions from TRISTAN were not in- 
cluded in the fitting procedure. There appears to 
be a discrepancy between the TOPAZ and AMY 
results. Indeed, this has been observed in a pre- 
vious study [f27f . 

4. RESULTS AND DISCUSSION 



Table 1 

HERWIG results (multiparton interactions on, 
not k t version). The x 2 per degree of freedom 
for a particular parameter set is given first for 
the combined fits (see text), and then for the sep- 
arate HERA and LEP fits. 



p.d.f. 




X 2 /D.o.F. 


HERA 


LEP 


WHIT 2 


1.6 


3.29 


4.24 


1.06 




1.8 


1.89 


2.27 


0.99 




2.0 


1.08 


1.11 


0.99 




2.2 


1.08 


1.07 


1.08 




2.4 


1.47 


1.60 


1.14 


GRV 92 


1.6 






1.03 




1.8 


1.47 


1.71 


0.89 




2.0 


1.35 


1.46 


1.09 




2.2 






1.17 


SaS Id 


1.4 


2.01 


2.30 


1.31 




1.6 






1.65 




1.8 


1.37 


1.16 


1.81 


SaS 2d 


1.6 


1.57 


1.81 


0.98 




1.8 


1.28 


1.40 


0.97 




2.0 






1.31 




2.2 


1.27 


1.25 


1.31 


LAC 1 


2.2 


8.70 


10.54 


4.35 




2.6 




2.14 





The results of the fits using HERWIG (not 
the kt version) are shown in table |l|, and those 
for PYTHIA in table |. Full details of all the 
fits, including the individual plots, can be found 



at http: / /www.hcp.ucl.ac.uk/^ jmb /HZTOOL / 
Figure |3j shows a summary of the results in graph- 
ical form. In most cases there is a clear, favoured 
value of p™ m , and overall - at least where reason- 
able fits arc found - there appears to be a favoured 



Table 2 

PYTHIA results (multiparton interactions on). 
The x 2 P er degree of freedom for a particular pa- 
rameter set is given first for the combined fits (see 
text), and then for the separate HERA and LEP 
fits. 



p.d.f. 




X 2 /D.o.F. 


HERA 


LEP 


WHIT 2 


2.4 


1.17 


0.91 


1.86 




2.6 


1.13 


0.99 


1.50 




2.8 


1.18 


1.07 


1.44 




3.0 


1.24 


1.12 


1.52 


GRV 92 


2.4 


1.21 


0.96 


1.89 




2.6 


1.24 


1.08 


1.66 




2.8 


1.24 


1.02 


1.72 


SaS Id 


1.4 


1.15 


1.15 


1.15 




1.6 


1.15 


1.19 


1.03 




1.8 


1.35 


1.38 


1.23 


SaS 2d 


1.4 


1.48 


1.32 


1.90 




1.6 


1.23 


1.21 


1.27 




1.8 


1.09 


1.08 


1.09 




2.0 


1.08 


1.07 


1.08 




2.2 


1.14 


1.11 


1.18 


LAC 1 


2.8 


6.18 


5.08 


10.0 



range of ~ 1.6 — 2.2 GeV. Note that although we 
do not constrain the p™ ln to have the same value 
for both LEP and HERA fits, when we combine to 
give an overall \ 2 we have taken those runs with 
the same p™ m . This is motivated by the idea that, 
as one moves toward low transverse momentum, 
two effects can render perturbative QCD inap- 
plicable. Not only does a s become large, but x 
becomes small, potentially leading to large ln(x) 
corrections. Our conjecture is that if a s effects 
are independent of the centre of mass energy, then 
p™ ln can be considered to be a universal param- 
eter. The low x effects are modelled by multiple 
interactions, and so the minijet model is "univer- 
sal", but different effects are seen depending on 
the cms energy and beam particle type. All the 
results shown have multiple interactions enabled. 
The agreement without multiple interactions is 
very poor. A particularly striking example of this 
is shown in figure ||, and it is low E^ et measure- 
ments such as this which are especially valuable 
for constraining multiple interaction models. 
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Figure 3. x 2 /DoF for combined LEP and HERA 
data samples (see text). 



The scaling of the overall normalisation gen- 
erally agrees between LEP and HERA for the 
good fits, and is in the range 1.2—1.3. In ad- 
dition, there is general consistency between dif- 
ferent datasets for both LEP and HERA. With 
regard to pdf sets, WHIT and GRV lead to bet- 
ter fits for HERWIG, whereas for PYTHIA there 
are no good fits to LEP data for sets other than 
the SaS sets. In no circumstances does it appear 
possible to obtain anything close to a good de- 
scription of the data with LAC1. 

The DIS tuned parameters for HERWIG do 
help, if only marginally, but there is no firm evi- 
dence yet that the HERWIG+fct modification im- 
proves matters. 

5. LINEAR COLLIDER 

The results of the tuning can be used to esti- 
mate photon-induced minijet backgrounds at fu- 
ture colliders. The best fits found using HERWIG 
have been extrapolated to linear collider ener- 



Figure 4. An example of the effect of multiple in- 
teractions. The data shown is ZEUS 1993 dijets 
data [2], along with the HERWIG predictions, us- 
ing the GRV 92 photon and pf in = 2.0, with 
the JIMMY multiple interactions model either en- 
abled or not. 



gies (y/s = 500 GeV), with the TESLA beam- 
strahlung and bremsstrahlung spectra included 
p8| ] . The parameters thus used were the WHIT 2 
photon pdf and a p™ m of 2.0 GeV with multiple 
interactions enabled. The overall normalisation 
was scaled by a factor of 1.4. The estimate for 
the minijet transverse momentum cross-section is 
shown in figure ||. 

The results are consistent at the 50% level 
with independent studies using PYTHIA and the 
DELPHI Monte Carlo ||. They indicate that 
the backgrounds will not be a concern for detec- 
tor occupancy and dosage, but that minijets po- 
tentially present a very significant background to 
physics owing to the high px tail. 
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Jet Transverse Momentum Spectra 
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